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We propose that axion-like particles (ALPs) with a two-photon vertex, consistent with aU astro- 
physical and laboratory bounds, may lead to a detectable signature in the spectra of high-energy 
gamma ray sources. This occurs as a result of gamma rays being converted into ALPs in the mag- 
netic fields of efficient astrophysical accelerators according to the "Hillas criterion", such as jets of 
active galactic nuclei or hot spots of radio galaxies. The discovery of such an effect is possible by 
GLAST in the 1-100 GeV range and by ground based gamma ray telescopes in the TeV range. 



o 
o 

o 

Q 



Oh' 
I 



(N 
> 
m 
o 

cn 

o 
l> 
o 



X 



PACS numbers: 98.70.Rz, 14.80.Mz 

Introduction — The Peccei-Quinn (PQ) mechanism fl*! 
remains perhaps the most compelling explanation of the 
CP problem of QCD. A new chiral [/pq(1) symmetry 
that is spontaneously broken at some large energy scale, 
/a, would allow for the dynamical restoration of the CP 
symmetry in strong interactions. An inevitable conse- 
quence of this mechanism is the existence of axions, the 
Nambu-Goldstone bosons of U{1)pq [1]. One of the most 
important phenomenological properties of the hypotheti- 
cal axion is its two-photon vertex which allows for axion- 
photon conversions in the presence of external electric or 
magnetic fields Q through an interaction term 



where a is the axion field, F is the electromagnetic field- 
strength tensor, F its dual, E the electric field, and B 
the magnetic field. The axion-photon coupling strength 
is quantified by 
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in Eq. ^ would lead to similar effects). This idea is 
somewhat similar to that discussed in the recent Ref . , 
but with some important differences. In that paper the 
authors considered the ALP parameters needed to fit 
PVLAS data [9] (as in other recently proposed gamma 
ray signatures of ALPs and assumed that the con- 

version of photons above ~ 10 TeV into ALPs takes 
place in the turbulent component of the galactic mag- 
netic field. Here, in contrast, we discuss the case in 
which the photon- ALP conversion occurs near or within 
the gamma ray sources. Interestingly, we find that, if the 
gamma sources are (or are hosted in) efficient astrophys- 
ical accelerators according to the "Hillas criterion" 
significant conversion can occur in ALP models which 
are fully consistent with all laboratory and astrophysi- 
cal constraints. In fact, the mechanism discussed here 
may offer the most practical way to detect ALPs over a 
significant range of masses and couplings. 

Photon-ALP conversion in gamma ray sources — As 
a consequence of the interaction of Eq. ([T|), ALPs and 
photons oscillate into each other in the presence of an 
(2) external magnetic field. For a photon of energy E^, the 
probability of converting into an ALP can be written [T^ 



(1) 



where a is the fine-structure constant and ^ is a param- 
eter of 0(1) depending on the details of the electromag- 
netic and color anomalies of the axial current associated 
with the axion field. In particular, this coupling is used 
by the ADMX experiment to search for axion dark mat- 
ter [1] and by the CAST experiment to search for solar 
axions The Peccei-Quinn axion has the important 

feature that its mass rria and interaction strength are in- 
versely related to each other and are connected to the 
measured properties of pions. One may, however, con- 
ceive of a more general class of particles whose coupling 
and mass are unrelated to each other. Such states are 
known as axion-like particles (ALPs). ALPs may mani- 
fest themselves in the propagation of photons in magnetic 
fields, either in laboratory or astrophysical environments, 
and may have potentially interesting astrophysical and 
cosmological consequences . 

In this letter, we propose another way to potentially 
detect ALPs, namely through their distortion of the en- 
ergy spectra of high-energy gamma ray sources (we note 
however that a light scalar particle coupling to F^^F^^'^ 



Pose = sin^(26l)sin2 




(3) 



where s is the size of the domain and B is the mag- 
netic field component along the polarization vector of the 
photon, which is assumed to be approximately constant 
within that domain. We have also defined an effective 
mixing angle 6 and characteristic energy E via 



sin2(26') = 
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where the effective ALP mass squared is = |r7i^ — Wpjl, 

cjpi = Ana n^lvfif^ is the plasma frequency, me the elec- 
tron mass, and rie the electron density. For the following 
considerations, it is useful to introduce the dimensionless 
quantities: g\\ — 307/10^^^ GeV~^, Bq, = _B/Gauss, 
Spc = s/parsec, m^ov = m/ixeSf, £goV = £/GeV. Re- 
cent results from the CAST experiment provide a di- 
rect bound on the ALP-photon coupling of gu < 8.8 for 
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iTT-a ^ 0.02 eV, nominally below the long-standing globu- 
lar cluster limit [7]. Note that 



Wpi = 0.37 X 10"'^ ^leV^/ne/cm~^ . 



(5) 



which means that in the interstellar medium (ISM) of 
the Milky Way, where Ue ~ 0.1 cm^'^, the effective mass 
of the ALP will not be smaller than m^cv ^ 10^^, in- 
dependently of how small nia is. For ultra-light ALPs 
{nia ^ 10~^^eV), the absence of gamma rays from 
SN 1987A yields a stringent limit oi gn ^ 1 [13] or even 
.911 < 0.3 [14]. But for 10"" eV < ma < lO'^eV, the 
CAST bound is the most general and stringent (bounds 
from ADMX, although stronger for some masses, assume 
that the axion is the Galactic dark matter). 

General properties — From Eqs. (j3l4p it follows that: 
(i) At energies below £, the mixing is small and a fortiori 
the conversion probability is small. Above this critical 
energy, the mixing is large, and a significant depletion 
probability might arise. In suitable units. 



^GcV = 



VcV 



OAgnBc 



(6) 



As we shall argue, when plugging in the previous formula 
typical astrophysical and ALP parameters, this critical 
energy naturally falls in the gamma ray energy range. 
Thus, the physics of light and weakly coupled ALPs nat- 
urally points to "f-rays as the most promising tool for 
discovery, (ii) A significant conversion into axions also 
requires that the argument of the oscillatory function in 
Eq. ([3]) is not too small, i.e. 



gajBs/2>l, i.e. 15 gn Spc ^ 1 • 



(7) 



The condition in Eq. ([7]) depends on the product B s, 
which also determines the maximum energy i?max to 
which sources can confine and thus accelerate ultra-high 
energy cosmic rays (UHECRs). This is known as the 
Hillas criterion [111 ], and for protons it writes 



9.3 X 102"cVBg spc 



(8) 



This connection between cosmic ray acceleration and 
Eq. ([7]) is important. Since UHECRs with energies of 
a few times lO^'' eV have been observed, environments 
where Bq Spc ^ 0.3 must exist in nature. This implies 
that couplings as small as gn c^i 0.2 might be probed, 
almost two orders of magnitude below present bounds. 

In the following, we shall examine in greater detail 
what are the signatures expected in 7-rays due to the 
ALP-photon conversion mechanism, what are the most 
promising sources to look at, and the perspectives for 
current instruments to probe the ALP parameter space. 

Spectral signatures — The qualitative signatures of the 
scenario considered here are remarkably robust, although 
the quantitative aspects are model dependent. The rea- 
son for this is twofold: (i) concerning particle physics, 
we ignore the fundamental mass and coupling parame- 
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FIG. 1: A typical power-law 7-ray spectrum (solid line) and 
its distortion for photon- ALP conversion with A = 1/3 and 
critical energies £ — 500 GeV (dashed lines) and £ = 2.5 TeV 
(dashed-dotted line). See text for details. 



3-D field configurations typically present in astrophysi- 
cal environments do not allow to apply naively Eq. Q 
for detailed quantitative predictions of the magnitude of 
the depletion. Nonetheless, the feature one should look 
for can be robustly parameterized as modification of the 
undistorted spectrum _Fo(i?-y) into a modified spectrum 
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1 + {£/E^r 



(9) 



ters nia and ga-^ 



the complicated (and unknown) 



where the constant A can be obtained in an idealized 
case from the oscillatory sinus function in Eq. As an 
estimate, an efficient conversion of an unpolarized pho- 
ton source resulting in a complete depolarization of the 
photon-ALP system would cause an average depletion 
of 33% of the initial photon flux (i.e. A ~ 1/3) above 
an energy £ determined roughly by Eq. ([6]). In Fig. [T] 
we plot the spectral distortion for ^ = 1/3 for the two 
cases £ = 500 GeV and £ = 2.5 TeV, and for a repre- 
sentative high energy gamma ray source. In particular, 
we assumed Fo(i?^) = k x , with r(= —2.2) and 
k consistent around 100 GeV with the flux of the blazar 
object Mkn 421 as reported by the MAGIC collaboration 
(15|. This object is a powerful emitter also observed by 
EGRET [11] at GeV energy. Note that to detect the sig- 
nature expected, one needs a wide dynamical range and 
sufficient statistics to detect a normalization shift of the 
typical power-law spectra at the level of 10-20%. Overall 
systematic errors in the energy scale, in the aperture and 
in the exposure are irrelevant to such a detection. Also, in 
variable sources the ALP signature should be impressed 
over the variable spectrum at all times, since it depends 
on the propagation and not on the emission of the pho- 
tons.) EGRET statistics and energy resolution are not 
sufficient to look for such a signature, and only recently 
ground-based gamma-ray telescopes have reached com- 
parable performances for powerful emitters within rea- 
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FIG. 2: Hillas diagram showing size and magnetic field 
strengths of astrophysical objects required to accelerate ultra- 
high energy cosmic rays (figure from Ref. [itI ] with permis- 
sion). The Hillas condition is closely related to the condi- 
tion for the efficient conversion of gamma rays into ALPs [see 
Eq. ©]. 



sonable exposure times. So, it would not be a surprise if 
such a signature had escaped detection so far, but would 
show up in the coming years thanks to the GLAST satel- 
hte detector and present and planned ground telescopes. 
In principle, if the astrophysical parameters were known, 
the amplitude of the depletion could be used to constrain 
Qa-y, see Eqs. (I3l7p . while the energy at which the effect 
is observed could be used to infer m, see Eq. ([6]). If only 
upper limits were available for B s, a lower limit on g^^ 
could be obtained, at least. Another prediction is that if 
a hint for an ALP would show up in a source at energy 
£i, then gamma emitters sitting in regions with similar 
values of B s should show a feature of similar amplitude 
at a characteristic energy related to £i only by the value 
of their field strength, see Eq. ([6]). 

Promising sources — To be consistent with existing 
bounds 6], an ALP should have a coupling gn ^ 9, which 
implies that Bq Spc ^ 10^^ must hold at the source. In 
Fig. [21 the Hillas diagram is shown, reporting the typi- 
cal B and s values for UHECR candidate sources. It is 
clear that virtually all the objects proposed as UHECR 
accelerators, from gamma-ray bursts to clusters of galax- 
ies, appear suitable for the search of ALP signatures. 
This is fortunate, since many observed (e.g. blazars) or 
expected (e.g. galaxy clusters) gamma ray sources are 
hosted in or near putative UHECR accelerators. How- 
ever, at least the compact sources on the Hillas plot are 
not likely the best candidates due to their higher densi- 
ties. Even a density of 10~^ g/cm'^, very low for terres- 
trial standards, would imply m ^ ImeV [see Eq. ([5])] 
and thus the condition E > £ cannot be satisfied in the 
energy range probed by gamma ray astronomy. This is 
a general, although qualitative, argument disfavoring too 
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compact (and presumably dense) sources as possible sites 
to observe photon- ALP mixing. Very promising sources 
are instead AGN jets and hot spots in radio galaxies such 
as Cygnus A and M87. For example, typical parameters 
for the hot spots of Cygnus A are Bq ~ 0.15 x 10^^, 
2 X 10'^ 18 1 , and similar numbers apply to the hot 
19[, which has been detected in the TeV 
In these environments, the quantity of Eq. ([7]) is 
near unity for gn ^ 0.3, while for propagation in our 
galaxy the same coupling would fail to satisfy that con- 
dition by more than one order of magnitude. 

A remark is in order. If it were proved that conser- 
vative estimates for the product B s of a detected 7- 
ray emitter satisfy Eq. ([7]) for gu ^ 9, then gamma 
observations would turn into powerful probes of ALP 
physics. But vice versa is not necessarily true: indeed, 
although some fits assuming synchrotron-self-Compton 
models seem to indicate that many detected gamma ray 
sources (see e.g. Refs. [H, HqI) fall short of the require- 
ment of Eq. ([7]) by one order of magnitude or more, it is 
important to remember that the ALP conversion feature 
depends on the properties of the environment crossed, 
not of the emitting region. In these cases, although a 
negative result can not be used to put significant bounds, 
a serendipitous discovery is by no means excluded. 

Exploring the ALP parameter space — One can easily 
estimate the range of ALP parameters observationally 
accessible. As we argued earlier, from the highest en- 
ergy UHECR observed the Hillas criterion suggests con- 
servatively that sites where Bq Spc ^ 0.3 must exist in 
nature. Once plugged into Eq. ([7]) (assuming equality), 
this implies that at least couplings as small as gu ~ 0.2 
may produce significant depletions in gamma-ray spec- 
tra. To deduce the range of masses which can be probed, 
we proceed as follows: (i) we neglect too compact objects 
on the Hillas plot and restrict our attention to the most 
promising range of astrophysical source sizes previously 
discussed, Spc ^ 10"** to 10^, deducing the corresponding 
field strength (see Fig. ^ ; (ii) we plug these values in Eq. 

Un- 
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(0) thus obtaining £qcy ^ 4 x 10" 
less TO ^ 0.1/xeV and the region is very compact (which 
is disfavored, as previously discussed), the transition en- 
ergy is expected to fall in the gamma ray band, confirm- 
ing again our initial, general remarks. Considering the 
energy range most interesting for GLAST (sub-GeV to a 
few tens of GeV) , we deduce that ALPs in the mass range 
W/icV ~ few X 10~^ 10^ can be probed. For ground 
based gamma ray telescopes such as HESS, MAGIC and 
VERITAS, which are sensitive to gamma rays in the ap- 
proximate range of 10^ to 10^ GeV, the photon- ALP con- 
version can also be significant for masses in the range 
of TO^oV ~ few X 10~^ to lO'^. Globally, we estimate 
the approximate range of parameters which could lead 
to observable effects as the region schematically shown 
in Fig. [3l along with the range excluded by CAST and 
the band preferred by QCD axion models The par- 
ticularly interesting region which overlaps with the QCD 
models band corresponds to 10 TeV transition energies 
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FIG. 3: The approximate range of ALP parameters which 
could lead to observable effects in gamma ray telescopes. Also 
shown are the parameters excluded by CAST and the band 
preferred for QCD axion models. 



with typical parameters of AGN cores. 

Summary — Space and ground-based gamma ray tele- 
scopes may have a chance to observe the effects of ALP- 
like particles (ALPs) through the mechanism of photon- 
ALP mixing. This mechanism leads to the depletion of 
gamma rays at high energies, resulting in a peculiar sig- 
nature in the spectra of gamma ray sources such as the 



jets of active galactic nuclei, the hot spots of radio galax- 
ies, or clusters of galaxies. If the astrophysical parame- 
ters at the sources were known sufficiently well, the mass 
and coupling of the ALP could be reconstructed or, in 
case of negative outcome, excluded. Without a detailed 
knowledge of the field strength and geometry around the 
gamma source a detection is still possible, but a neg- 
ative result cannot easily be translated into interesting 
exclusion plots in the ALP parameter space. We believe 
that the mechanism discussed here should be thought of 
as an example of the opportunities that the new gen- 
eration of gamma ray telescopes will offer for studying 
fundamental physics. Note that gamma ray telescopes 
may observe the effects of ALPs with parameters which 
are exceedingly difficult to explore otherwise, as shown 
in Fig. [3) An appearance experiment such as CAST is 
actually sensitive to g'^^, and thus would require an im- 
provement of five to six orders of magnitude in sensitivity 
to cover the entire range of parameters considered here. 
However, a detection of the kind described here could be 
confirmed, in part of the parameter space, by other astro- 
physical techniques, such as that suggested in Ref. [2l| . 
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